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Abstract. A new damage model based on a micromechanical analysis of cracked [±#°/90°] s 
laminates subjected to multiaxial loads is proposed. The model predicts the onset and accumu- 
lation of transverse matrix cracks in uniformly stressed laminates, the effect of matrix cracks 
on the stiffness of the laminate, as well as the ultimate failure of the laminate. The model also 
accounts for the effect of the ply thickness on the ply strength. Predictions relating the elastic 
properties of several laminates and multiaxial loads are presented. 
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1 INTRODUCTION 

The aerospace industry is committed to improve the performance of aircraft whilst reducing 
emissions and weight. Such a goal can be achieved by the use of advanced polymer-based com- 
posite materials, that have excellent properties for aerospace applications, such as low density, 
and fatigue and corrosion resistance. 

The design procedure used for advanced composite structures relies on a ’building-block’ 
approach [1 J, where a large number of experimental tests are performed throughout the product 
development process. The use of improved analytical or numerical models in the prediction 
of the mechanical behavior of composite structures can significantly reduce the cost of such 
structures. Such models should predict the onset of material degradation, the effect of the non- 
critical damage mechanisms on the stiffness of the laminate, and ultimate structural failure. 

Strength-based failure criteria are commonly used to predict failure in composite materi- 
als L2J-L7J. Failure criteria predict the onset of the several damage mechanisms occurring in 
composites and, depending on the laminate, geometry and loading conditions, may also predict 
structural collapse. 

In multidirectional composite laminates, damage accumulates during the loading process. 
Final failure occurs as a result of damage accumulation and stress re-distribution. The ulti- 
mate failure load is higher than the initial damage predicted by strength-based failure criteria. 
Furthermore, stress- or strain-based failure criteria cannot represent size effects that occur in 
quasi-brittle materials [8J. 

Simplified models, such as the ply discount method where some scalar components of the 
stiffness tensor are reduced to approximately zero when damage is predicted, are often used by 
the industry to predict ultimate laminate failure. However, these methods cannot represent with 
satisfactory accuracy the progressive reduction of the stiffness of a laminate as a result of the 
accumulation of matrix cracks. 

Constitutive laws based on Continuum Damage Mechanics (CDM) have been proposed to 
predict the material response, from the onset of damage up to final collapse [9]-[14J. Although 
the existing CDM models can predict accurately the evolution of damage, they usually rely 
on fitting parameters that need to be measured at laminate level, such as the critical values of 
thermodynamic forces [14J. 

Three-dimensional models based on CDM use material properties measured at the ply level. 
Damage mechanisms such as transverse matrix cracks are represented by strain- softening con- 
stitutive models and bands of localized deformation. The implementation of strain-softening 
models in finite element models causes convergence difficulties during the iterative solution 
procedure. Furthermore, strain-softening models require regularization techniques to provide 
mesh-independent solutions [15, 16J. 

Alternative methods based on the combination of elastic analysis of cracked plies and finite 
Fracture Mechanics provide the basis for an accurate representation of the response of compos- 
ite materials [ 1 7 J-[28J . Micromechanical models have been developed to predict the initiation 
and evolution of transverse matrix cracks under either in-plane shear or transverse tension. Gen- 
eralizations of these models are required for the more usual case of multiaxial loading. 

In order to predict ultimate failure, a micromechanical model needs to be used in combination 
with a fiber failure criterion. Furthermore, a general methodology to predict laminate strength 
should be able to predict matrix cracking under high values of transverse compression. These 
damage mechanisms usually correspond to the final failure of laminates under uniform states of 
stress [6J. 
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The objective of this work is to define a new damage model based on micromechanical 
models of transverse matrix cracks to predict the onset and evolution of transverse matrix cracks 
under multiaxial loading. A new constitutive model is derived based on the thermodynamics 
of irreversible processes. The model proposed herein represents transverse matrix cracks as 
distributed damage. The model is able to predict the onset and propagation of matrix transverse 
cracks under multiaxial loading as well as the final failure of uniformly stressed laminates where 
a periodic distribution of transverse matrix cracks can be assumed. 

2 MICROMECHANICAL MODEL 

The proposed continuum damage model is based on two major components: a set of stress 
based failure criteria and a micromechanical model of transverse matrix cracking in multidirec- 
tional laminates. 

The failure criteria define the onset of transverse matrix cracking, i.e. the activation of the 
damage variables. A micromechanical model of transverse matrix cracks is required to define 
the evolution of the damage variables. 

Several micromechanical models of transverse matrix cracks that have been proposed in the 
literature [19, 28 J can be used within the framework developed here. The micromechanical 
model proposed by Tan and Nuismer [17, 18J accounts for the effects of the adjoining plies on 
the homogenized elastic properties of a cracked ply. This micromechanical model is the basis 
of the developments presented in this work. 

Using the assumption of generalized plane strain, Tan and Nuismer [17, 18J developed a 
model able to relate the density of transverse matrix cracks in a central 90° ply to the homoge- 
nized elastic properties of that ply. The model developed by Tan and Nuismer was used for the 
prediction of the evolution of transverse matrix cracks under either in-plane shear or transverse 
tensile stresses. 

The laminates under investigation are symmetric and balanced with a [±0/9O°] s layup con- 
taining a periodic distribution of transverse matrix cracks, as shown in Figure 1 . The micro- 
mechanical analysis of a balanced symmetrical laminate requires the division of the laminate 
in two sub-laminates: the 90° layers in the middle layer (sublaminate 1), and the outer plies 
(sublaminate 2). 



Figure 1 : One-quarter of the laminate. 


Sublaminate 1 is taken as transversely isotropic. The outer layers, sublaminate 2, may con- 
tain several layers with different fiber orientations but must be orthotropic in the laminate global 
system (x, y and z coordinates). 
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2.1 CONSTITUTIVE TENSOR 

The stiffness tensor of the balanced and symmetric laminate shown in Figure 1 can be written 
as a function of the density of transverse matrix cracks (1/L) in the sublaminate 1: 


ML) 


A{L) ii A(L) 12 0 

A(L) 21 A(L ) 2 2 _ o 
0 0 A(L)qg 


( 1 ) 


where the terms Aij(L),i,j = 1, 2, 6 are obtained from the Tan and Nuismer model [18J. 

The undamaged stiffness matrix Q of the laminate is obtained from lamination theory using 
the undamaged stiffness matrix of sublaminate 1, Q n) , and the stiffness matrix of sublaminate 

2,QW: 


<?«’ = ^ (s<3« - t^of) (2) 

with h = + h^ 2 \ 

Assuming that the degradation due to the transverse matrix cracks only occur in sublaminate 
1, the damaged stiffness tensor of laminate 1 is given as: 

^ - h^Qf) (3) 

where (L) (i, j = 1,2,6) are the scalar components of the stiffness tensor of sublaminate 1. 
These components are a function of the distance between transverse matrix cracks (L. defined 
in Figure 1). 

Having defined it is possible to calculate the effective transverse modulus, Poisson 

ratio, and shear modulus of the 90° ply: 


E^(L) = E<i\L) = 4;>(L) - [ ( J 

A 22 \L) 

(4) 

Ml \L) = My (L) = ^ L) 
’ a£\l) 

(5) 

G { S(L) = G«(L) = A&\L) 

(6) 


The quotient ^ is not a function of damage. This observation is in agreement with several 
other models, such as the ones proposed by Laws et al. [23, 24J, and Nguyen [25]. 

The plane stress compliance tensor of the damaged sublaminate 1, H (1 ), only contains two 
components that depend on the density of transverse matrix cracks: H M (L) and Hqq (L). The 
tensor H (1 ), is established as a function of the density of transverse matrix cracks, L, as: 


H (1 ) = 


j_ 
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with: 


«£ (l) = 


E. 


(i) 


A’ 1 1 1 

W A[\\l)a$(l) - \a^(l) 


h£ (l) = 


G { S(L) 


A m 


*66 


(L) 


( 8 ) 

(9) 


The evolution of the homogenized elastic properties of the cracked ply (sublaminate 1) as 
a function of (3 = 1/(2 L) is shown in Figure 2 for a [±25°/903] s laminate with the material 
properties shown in Table 1 , 


Ei (GPa) E 2 (GPa) 

Gi 2 (GPa) 

v l2 

163.4 11.9 

6.2 

0.30 


Table 1: Elastic properties of typical carbon/epoxy composite. 



Figure 2: Effective elastic properties of the 90° ply as a function of the density of matrix cracks. 

Figure 2 shows that the Young modulus of sublaminate 1, E A , is not affected by the presence 
of transverse matrix cracks. Furthermore, the curves E 2 (/3)/E 2 and v 2 i(/3)/v 2 \ are coincident 
(Figure 2). 

2.2 ONSET OF MATRIX TRANSVERSE CRACKS 

The onset of transverse matrix cracking in a ply under the combined effect of in-plane shear 
stresses and transverse tensile stresses needs to be predicted using an appropriate failure crite- 
rion. The failure criterion should be established in terms of the actual strengths of a ply when it 
is embedded in a multidirectional laminate (in-situ strengths). 
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2.2.1 In-situ strengths 

A methodology to predict the onset of matrix transverse cracking must be able to predict 
the in-situ strengths of a ply. Both the transverse tensile and in-plane shear strengths of a ply 
embedded in a multidirectional laminate are higher than the corresponding values measured in 
unidirectional laminates [29J. 

The thick models described in [30J are used. The tensile transverse in-situ strengths of 
sublaminate (1) is [30J: 


For a thin embedded ply: Yj = 



For a thick ply: Y T = 1.12\/2Yj d 


( 10 ) 

(ID 


where Y^ d is the tensile transverse strength measured in a unidirectional test specimen, t is the 
thickness of sublaminate 1, G\ c is the mode I intralaminar fracture toughness and A 22 is defined 
as: 


v. 


A° 22 = 2 ( ^ 


21 


E 2 E\ 

The in-situ shear strengths are obtained as [29J: 


( 12 ) 


St = 


(1 + y0Gf 2 ) ^ 


(13) 


3xG*i2 

where x is the shear response factor defined in [29 J, and the parameter </> is calculated according 
to the configuration of the ply: 


For a thick ply: 
For a thin ply: 


12 (jgf 

G,2 

48G'n c 

Tit 


udl ^ 


+ 18 * (S? ) 


(14) 


where S^ d is the shear strength measured using an unidirectional test specimen and 6' llc is the 
mode II fracture toughness. 


2.2.2 Failure criterion for the prediction of transverse cracking under multiaxial loading 

In general, a ply represented by sublaminate 1 in Figure 1 is subjected to transverse tensile 
stresses and in-plane shear stresses. Under pure in-plane shear or pure transverse tension, the 
onset of transverse matrix cracking is predicted by comparing the components of the stress 
tensor with the respective in-situ strengths (defined in the previous section). 

Under multiaxial loading, it is necessary to use a failure criterion to predict the onset of 
matrix cracking. The LaRC04 [4, 5] failure criteria are a function of the components of the 
stress tensor and in-situ strengths. For transverse tension, the criterion used is: 

(1 - 0)?r + 3 (t?) 2 + (-?) 2 - 1 < 0 with ci 22 > 0 (15) 

it r t ol 
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where g — Gn c is the mode II component of the fracture toughness associated with matrix 
transverse cracking. 

Under moderate values of transverse compression, the cracks propagate along the laminate 
thickness direction, and the LaRC04 failure criterion is: 


7T (1 012 1 + V L(J 22) - 1 < 0 with a 2 2 < 0 
where // L is the coefficient of longitudinal influence defined in [4, 5J. 


(16) 


2.3 Evolution of matrix transverse cracks under multiaxial loading 
2.3.1 Transverse tension 

Tan and Nuismer [18J obtained a closed-form expression that defines the evolution of trans- 
verse matrix cracks under uniaxial stress states (either transverse tension or in-plane shear). 

To predict failure under multiaxial loading, i.e. when the lamina is simultaneously under 
tensile and in-plane shear strains, e xx and 7 xy respectively, it is necessary to derive a relation 
between the density of transverse matrix cracks and the applied multiaxial strain state. It is 
assumed that the relation between the tensile and shear strains is constant throughout the loading 
history: 


k = with e xx > 0 

&XX 

where k is the multiaxial strain ratio. 

Consider a ply with crack density L, as shown in Figure 3 a). 


(17) 









< > 



a) Crack density 2L 



Figure 3: Progression of transverse matrix cracks. 


The strain energy in a laminate cell of length 2 L subjected to transverse tension and in-plane 
shear just prior to fracture, U 2 L, can be established as a function of the strain tensor and of the 
crack density as: 


U 2L = 2 hLb [E,(L)el x + A^L)^} (18) 

where E X (L ) is the axial stiffness of the laminate with a crack density of 1 / (2L) and b is the 
specimen width; A ee (L) is the laminate effective shear stiffness corresponding to a crack density 
of 1/(2L), 2 h and 2 L are the laminate thickness and the distance between two consecutive 
transverse matrix cracks, respectively. 

After the propagation of transverse matrix cracks, Figure 3 b), the strain energy in the original 
unit cell of length 2 L is: 


U L = 2 hLb 




+ A 


66 



(19) 
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where E x (L/2) and A m (L/2) are respectively the axial stiffness and the laminate effective 
shear stiffness corresponding to a crack density defined by L. 

The energy required to generate a new matrix crack in a ply equals the loss of strain energy 
of the laminate [18J. Therefore, the difference between equation (18) and equation (19) is equal 
to the energy released by the sublaminate 1 : 


AU 


U 2L -U L = 


2 hLb 


E X (L) — E x 



E 2 + 

^XX 1 


Am(L) — A 66 




( 20 ) 


AU = 2 h (1) bG c = 2 h w bGj + 2 h w bG n (21) 

where G c is the mixed-mode fracture toughness of sublaminate 1 under tensile (mode I) and 
shear (mode II) loading. From (20) and (21): 


E X (L) — E x 



£ 2 ~\~ 
^XX ' 


^66 (L) — A 66 




hWGc 

hL 


( 22 ) 


Using the definition of k given in (17), equation (22) can be re-written as a function of the 
strains: 


or: 


{(^(L) -£*(§)) 


+ Kf 


— A 66 



h^G c 

hL 



E X (L ) — E x 



+ 


^66 {L) — A m 




h&G c 

hL 


(23) 


(24) 


The relations between the normal and shear strains and the crack density are obtained as: 


IhWGr. 


&XX 


hL [E X (L) - E x (f )] + k 2 [A m (L) - A 66 (f)] 


(25) 


Ixy 


<hWG c 

hL \E X (L) — E x + h 2 [Aq 6 (L) — A ee (|f)] 


KjEst 


(26) 


Equations (25 ) and (26), are established as a function of the mixed-mode fracture toughness 
G c that needs to be defined. 

The mixed-mode fracture toughness is defined in terms of the mode I and mode II compo- 
nents as: 


From (27 ) and (22 ): 


G c — Gj + Gjj 


h^Gj 

hL 

hP>G n 

hL 


E X (L) — E x 



Am{L) — Aq & 




(27) 

(28) 
(29) 


8 


Pedro P. Camanho, Joan A. Mayugo, Pere Mailin' and Carlos G. Davila 


Dividing equations (28) and (29) by (23) and (24), respectively: 


Gj _ [- E X (L ) - E x (f )] 

" G c [£ X (L) - E x (f )] + k 2 [A 66 (L) - A 66 (§)] 

G// _ ft 2 [^66 (■£') — 2L66 (|f)] 

= |X(L) - E x (f )] + « 2 \A m (L) - A 66 (f )] 


A/ and A// are the ratios between the mode I and mode II components of the energy release 
rate and the mixed-mode fracture toughness. 

The criterion proposed by Hahn [31 ] for the prediction of transverse matrix cracking under 
transverse tensile and in-plane shear loads is used: 


(1 ~9) 



. G l 
+ 9yr~ 

Oric 


+ 


G\\ 

Guc 


Substituting equations (30) and (31) into (32) gives: 


1 


(32) 


(1 ~9) 


ArGr 


G 


+ 9- 


ArG, 


Ic 


G 


c + { — 


A t )G c 


Ic 


G 


= 1 


lie 


The positive real solution of (33) is: 


(33) 


G c — G Uc + 


A i (G*ic — G\\ c y 
2 Gl r 


'1 + 


G, r G 


Ic Vi He 


A I (Gi c — Gr'nc)" 


(34) 


where Aj depends on the density of transverse matrix cracks, 3. and on the multiaxial strain 
ratio, k (equation (30)). 

Figure 4 shows the relation between the crack density 3 and s xx for different multiaxial 
strain ratios. A [ ±25° / 902] , s laminate with the elastic properties shown in Table 1 is used in the 
predictions. 


8 XX 



Figure 4: Relation between applied strain ( e xx and k) and density of matrix cracks. 
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The effects of multiaxial strain states on the crack density are clearly shown in Figure 4, 
It can be observed in Figure 4 that the density of transverse matrix cracks increases with the 
multiaxial strain ratio for a fixed value of e xx . 

2.3.2 Transverse compression 

Transverse matrix cracks created by a combination of transverse compression and in-plane 
shear close under the effect of the compressive transverse stress. When a crack closes, its faces 
can transmit normal tractions but shear tractions may cause slippage between the crack faces. 
Therefore, it can be assumed that transverse matrix cracks only affect the shear stiffness of a 
laminate. 

Following the procedure described in the previous section, the energy released by the sub- 
laminate 1 is: 


AU 


2hLtrfl y 


Am(L) — Aq 6 



2h {1) bG llc 


The relation between the shear strain and the crack density is obtained as: 


(35) 


hL [Aqq(L) — Aqq (|t)] 

3 MICROMECHANICS-BASED DAMAGE MODEL 
3.1 Constitutive model 

A damage model able to represent the onset and accumulation of a periodic distribution of 
transverse matrix cracks should yield a compliance tensor similar to the one obtained from 
the micromechanical model, equation (7 ). An appropriate damage model can be developed by 
defining the Gibbs free energy per unit volume, To, as: 



= - 


a 


E i ^ 


a, 


22 


+ 


cr 


12 


^21 ^12 
E‘2 Ei 


(1 — d 2 ) E 2 (1 — d 6 ) G 12 
+ (a n<Tn + a 22 (J 22 ) AT + + (3 22 <j 22 )AM 


&11&22 


+ 


(37) 


where E u E 2 , v 12 , v 2 \ and G V2 are the in-plane elastic orthotropic properties of a unidirectional 
lamina. The subscript 1 denotes the longitudinal (fiber) direction and 2 denotes the transverse 
(matrix) direction. d 2 and <7 fj are damage variables associated with transverse matrix crack- 
ing. a n and a 22 are the coefficients of thermal expansion in the longitudinal and transverse 
directions, respectively. /3 U and 8 22 are the coefficients of hygroscopic expansion in the lon- 
gitudinal and transverse directions, respectively. AT and AM are respectively the changes in 
temperature and moisture content from the stress-free state. 

The proposed model distinguishes between active (d 2+ ) and passive damage (d 2 -) variables, 
corresponding to the opening or closure of transverse matrix cracks under load reversal respec- 
tively. To account for the active damage under either opening or closure, the following equation 
is proposed: 


d 2 


d 2 + 


(cr 22 ) 

I <^22 1 


+ d 2 — 


{-(y 22 ) 

I <^22 1 


(38) 


10 


Pedro P. Camanho, Joan A. Mayugo, Pere Maimf and Carlos G. Davila 


where the McCauley operator (x) is defined as (x) := \ (x + \x\). 

The damage variables d 2 and d 6 can be related to the density of transverse matrix cracks 
using the micromechanical model previously described. 

The constitutive model is obtained as: 


£ =-xr^~ — H : a+ATa + AM/3 (39) 

da 

where the compliance tensor H is: 


H 


da 2 


1 V21 

E\ E 2 

V12 1 

Ei (l—d2)E2 
0 0 


0 

0 

1 

(l—de)Gi2 


(40) 


The compliance tensor is established in terms of the damage variables and is similar to the 
compliance tensor derived in the micromechanical model. 

Using (7) and (40) the damage variables d 2 and d 6 can be expressed in terms of the crack 
density (3 as: 


d 2 + — 1 — 


1 

E 2 H 22 (/3) 


(41) 


da — 1 


1 

Gi 2 H 66 {(3) 


(42) 


The functions H 22 (/3) and 77 66 (/3) are obtained from equations (8) and (9) with f3 = 1/(2L). 
The condition of positive energy dissipation is established as a function of the thermody- 
namic forces, Y, and as a function of the time derivative of the damage variables, d, and it is 
given by: 


= Yd = 


d^c 

dd 


• d > 0 


a. 


22 


2(1 — G ? 2 +) ^2 


d'2+ + 


a 


12 


2(1 — do) G 12 


-d 6 >0 


(43) 


where S is the rate of energy dissipation per unit volume. 

The condition of positive dissipation, when interpreted from a micromechanical point of 
view, establishes that the density of transverse matrix cracks can only increase or remain con- 
stant. 


3.1.1 Damage activation functions 

Transverse matrix cracks are predicted using two scalar functions, F k (k = 2+, 2—), estab- 
lished in terms of the effective stress tensor d t , and of the damage threshold value, r tm . 

F k ■= 4>k (d f ) - r* < 0 (44) 

where t is the current time. The initial threshold value, r°, is equal to 1, and the following 
condition must be satisfied in order to fulfill the Second Principle of Thermodynamics: 

f > 0 (45) 
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Damage onset occurs when any of the functions 4>k (cr*) reaches the initial damage threshold 
value of 1. The functions (pk (<7*) used to predict transverse matrix cracking are based on the 
LaRC04 failure criteria previously proposed by the authors [4, 5J. The damage activation func- 
tions are used to predict the onset of transverse matrix cracks lying in ply thickness direction, 
as shown in Figure 1 . 

Transverse matrix cracks lying in the direction of the ply thickness occur under transverse 
tension (022 > 0), or under moderate values of transverse compression and in-plane shear 

(<722 < 0 ). 

For high values of transverse compression, the matrix fractures lie along a plane that is 
inclined at an angle a to the direction of the ply thickness. Inclined fracture planes caused 
by transverse compression induce delamination between the plies. This damage mechanism is 
usually catastrophic in uniformly stressed composites where local redistribution to more lightly 
loaded regions of the structure cannot occur [6J. Therefore, laminate failure can be assumed to 
occur when matrix cracking under high values of transverse compression is predicted. 

The damage activation function used to predict matrix cracking under transverse tension 
(a ! 22 > 0) and in-plane shear is defined as: 

F 2+ := 0 2 + (cf*) - r l < 0 (46) 

with: 


<4?) 

The damage activation function used to predict matrix cracking with a = 0° under moderate 
values of transverse compression (at 22 < 0) and in-plane shear is defined as: 

F 2 _ := 0 2 . (cr*) - r* < 0 (48) 

with: 

02- (a*) = (\d-[ 2 \ + ri L al 2 ) (49) 

a=0° Dl 

The equations proposed can predict accurately transverse matrix transverse cracking. The 
failure envelope, or initial elastic limit, predicted using equations (46) and (48), and the corre- 
sponding experimental data are shown in Figure 5 . 
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Figure 5: Comparison between predictions and experimental data 1 32 1 and definition of ply or laminate failure 
domains. 

Figure 5 also shows the domain of validity of the model: the damage model is only defined 
for fracture angles a = 0°. For a > 0°, the laminate is assumed to fail. 

3.1.2 Damage evolution functions 

Using the principle of strain equivalence, the effective stress tensor can be written as a func- 
tion of the strain and undamaged stiffness tensors: 

a = C°: e (50) 

The strain tensor is established in terms of the density of transverse cracks and the multiaxial 
strain ratio k, equations (25-26). Therefore, all terms in the damage threshold functions can be 
formulated as a function of the density of transverse matrix cracks and of the multiaxial strain 
ratio: 


F k = f(/3\ k*) 


(51) 


For a given state of strain at time t, the multiaxial strain ratio is a dependent variable that can 
be easily defined using equation ( 17) in material coordinates: 


K t-Jk 

K - .t 


£ 


11 


(52) 


The density of transverse matrix cracks is a state variable. Therefore, it is necessary to define 
an evolution law subjected to thermodynamic restrictions. 

The first condition to be satisfied is the requirement of positive dissipation: 


$ > 0 


(53) 
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Furthermore, it is necessary to define the conditions for the evolution of the elastic domain: 


F k {a\p) <0 


(54) 


f3F k (Z\p t ) = 0 (55) 

Equations (53 )-( 55 ) are the Kuhn-Tucker conditions which ensure a consistent evolution of 
damage during loading and load reversals. 

The evolution laws of the state variables are defined as: 


a 


r 


da . da 
de de 


dr ■ dr 

dp^ + d^ K 


de ■ de . 

dpP + dK K . 


C° : 



de 

+ d~. K 


(56) 

(57) 


For a given loading state, the damage consistency condition must be applied to define the 
evolution of the internal variables. The consistency condition is defined as: 


From (46) and (48): 
Using (59) in (58): 


F k — 0 => F k = 


OF, 


k ~ 


dd 


dF k 

a + -=—f = 0 
dr 


dF± 

dr 


(58) 

(59) 


dF k ■ dF k 

F k = -^-.a-r = -^:C°-.e-r = 0 
da da 


Using equations (56) and (57) in (60): 


dFi 


f de 


de 


F k = —F:C°:\ — /3 + — k - — -/? + — k =0 


da 


\d/3 dn 


( dr 


dr 


\d/3 dn 


(60) 


(61) 


Taking into account that the micromechanical model proposed assumes a constant loading 
ratio, k = 0: 


F k 


dFk f c o 
da \ ' df3 



(62) 


The density of of transverse matrix cracks, /3, is calculated from the integration of equation 
(62) using a numerical method, such as the return-mapping algorithm. 


4 LAMINATE FAILURE 

In the presence of stress concentrations the onset of fiber localized failure does not cause 
immediate structural collapse. Experimental results [33, 34J show that structural collapse is 
caused by the the progression of fiber fracture. Therefore, it is necessary to use a damage model 
that accounts for the stress re-distributions caused by fiber fractures. 

The model proposed here assumes that for uniformly stressed laminates the onset of fiber 
fracture and delamination caused by high compressive transverse stresses triggers structural 
collapse. Therefore, laminate final failure is predicted when fiber failure or matrix cracking 
with a 7^ 0° occurs. 
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4.1 Fiber failure 

The criterion for fiber fracture under longitudinal tension (an > 0) is defined as [4, 5 J : 

FIi + := ^ — 1 < 0 (63) 

Ax 

where X T is the ply tensile strength in the longitudinal direction. 

The LaRC03 failure criterion for fiber kinking under longitudinal compression (an < 0) is a 
function of the components of the stress tensor in a frame representing the fiber misalignment, 
<#° L4, 5J: 

a© = a n cos 2 p + a 2 2 sin 2 ip + 2 |ai 2 | sin p cos p 

a© = a n sin 2 + cr 22 cos 2 p — 2 |ai 2 1 sin p cos p (64) 

a 1 = — cTn sin p cos p + a 22 sin p cos p + |a i2 1 (cos 2 p — sin 2 p) 

where the fiber misalignment angle p is defined as [4J: 


<P 


Cl2 1 + {G\2 ~ Xq) p C 
G 12 + on — a 22 



(65) 

( 66 ) 


where Xq is the ply compressive strength in the longitudinal direction. 

Depending on the sign of the in-plane transverse stress \ the criteria for fiber kinking 
(an < 0) are: 


FIi 


or Fin 



1 < 0, 


+ 


( m ) ^ n 

0-22 < 0 


a 


(m)' 

12 


5 t . 


1 < 0, 


a. 


(m) 


22 — 


> 0 


(67) 

( 68 ) 


4.2 Matrix failure with a ^ 0° 

The failure criteria for matrix cracking under transverse compression (a 22 < 0) and in-plane 
shear and a ^ 0° are defined as [4, 5J: 


Fb - : =S) + (£) -^ 0 .Pn!< Yc (69) 

/ (m)T\ 2 /(m)L\ 2 

Fl2 ‘ ;= W + W " 1 -°’ kul>Yc <?0) 

where the effective shear stresses in the fracture plane are defined as: 
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rj — {\t t \+ rf o n cos 9) (71) 

T e = (\r L \ +rj L cr n sin 9) (72) 

where 6 = tan -1 ( ) . The components of the stress tensor on the fracture plane are given 

\ (T22 sin a J r r o 

by [4, 5 J : 

{ On = 022 COS 2 Ct 

r T = — cr 22 sin ct cos ct (73) 

t l = cr 12 cos a 

The terms t™ 7 and r"' 7 ' are calculated from equations (71 )-(73 ) using the relevant compo- 
nents of the stress tensor established in a frame representing the fiber misalignment, equation 
(64). The angle a is determined by maximizing the failure index FI 2 _ (69-70) using a simple 
iterative procedure. 

The coefficients of transverse and longitudinal influence, r/ T and // 7 , respectively, are [4, 5 J : 


tan 2ct 0 

L S L cos2a 0 

V = - v 2 

Yc cos^ a o 

with O'o rs 53°. Yq is the ply compressive strength in the transverse direction. 

5 EXAMPLES 

The present damage model can be used in combination with classical lamination theory 
using stand-alone codes. Alternatively, the damage model can be implemented as a constitutive 
subroutine for the finite element method. 

The damage model was implemented using a commercial symbolic computing software. The 
model was verified by calculating the response of several glass-epoxy laminates under uniaxial 
and multiaxial loads: [±45°/90^, [0^/90^] s , [0§/90§] s and [0^/90°],. 

In all calculations performed, the ply thickness was taken as 0.144mm, and the temperature 
difference from the stress free condition was -100°C. The coefficients of thermal expansion in 
the longitudinal and transverse directions are «n = 7.43 x 10” 6 /°C and a 2 2 = 22.4 x 10 _6 /°C, 
respectively. The remaining material properties used are shown in Tables 2 and 3 . 


(74) 

(75) 


E, (GPa) 

E 2 (GPa) 

G 12 (GPa) 

G 23 (GPa) 

Vl2 

V23 

X (10 -8 MPa~ a ) 

44.7 

12.8 

5.8 

4.5 

0.30 

0.42 

2.0 


Table 2: Elastic properties of glass-epoxy [ 20 1 
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Y" d (MPa) 

S; id (MPa) 

Gic (N/mm d ) 

Guc (N/mm ') 

40.0 

73.0 

0.20 

0.40 


Table 3: Strengths and fracture toughnesses of glass-epoxy [ 20 1 . 


The predicted response of [±40° /QO^ and [O^/OO^s, is shown in Figure 6 for two values of 
the multiaxial strain ratio: k = 0 and k = 10. 


Ex(P)/Ex 



Figure 6: Relation between laminate modulus and applied strain for [±45 o /90|] s and [0§ / 90^] s laminates. 


Figure 6 shows that the rate of degradation of the elastic properties of the laminate is higher 
when the axial stiffness of the outer sublaminate decreases. The effect of multiaxial loading 
is also clear in Figure 6: as expected, the application of shear strains leads to a reduction of 
the extension corresponding to the onset of transverse matrix cracks and to a higher rate of 
degradation of the elastic properties of the laminate. 

Figure 7 compares the response of [O^/OO^].;, [O^/ilO^s and [02/90°] s laminates for n = 0 
and k — 2. 


Ex(p) /Ex 



Figure 7: Relation between laminate modulus and applied strain for [0§ / 90^] [01 / 90§] s and [O 2 / 9 O 0 ] ,, laminates. 
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The in-situ effect is shown in Figure 7 : for k = 0, the strain corresponding to the onset 
of matrix cracking of the [02/90°] s laminate is 1.9 and 1.3 times higher than the strains of the 
IO2/9O4]* and [O2/9O2] s laminates respectively. Furthermore, the strain at the onset of matrix 
cracking of the [Oy/OO 0 ],, is 2.7 times higher than the ultimate transverse strain measured in an 
unidirectional test specimen. 

6 CONCLUSIONS 

A new, micromechanics-based, continuum damage model able to simulate the onset and 
propagation of transverse matrix cracks and final laminate failure is proposed. The model is 
applicable to [±#°/90°] s laminates, under multiaxial loading and uniform stresses or small 
stress gradients. 

The model uses ply properties and does not require any tests performed at the laminate 
level to identify damage onset and evolution functions. The onset of damage is predicted using 
failure criteria and damage evolution laws are established from the micromechanical analysis 
of cracked plies. 

The onset and accumulation of transverse matrix cracks are represented as a distributed dam- 
age mechanism. The onset of localization, which is triggered by either fiber fracture or matrix 
cracking with a 7^ 0, is assumed to cause a structural collapse. 

The predictions show that the rate of degradation of the elastic properties increases when 
the stiffness of the outer sublaminate decreases. Decreasing the thickness of the 90° plies also 
increases the degradation rate. 
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